XMM-Newton Observation of the Merging Galaxy Cluster A1644 
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ABSTRACT 

An XMM-Newton imaging spectroscopy analysis of the galaxy cluster A1644 is presented. 
A1644 is a complex merging system consisting of a main and a sub cluster. A trail of cool, 
metal-rich gas has been discovered close to the sub cluster. The combination of results from 
X-ray, optical, and radio data, and a comparison to a hydrodynamical simulation suggest that 
the sub cluster has passed by the main cluster off-axis and a fraction of its gas has been stripped 
off during this process. Furthermore, for this merging system, simple effects are illustrated which 
can affect the use of clusters as cosmological probes. Specifically, double clusters may affect 
estimates of the cluster number density when treated as a single system. Mergers, as well as cool 
cores, can alter the X-ray luminosity and temperature measured for clusters, causing these values 
to differ from those expected in equilibrium. 

Subject headings: cooling flows — galaxies: clusters: general — galaxies: clusters: individual (Abell 
1644) — intergalactic medium — X-rays: galaxies: clusters — cosmology: observations 



1. Introduction 

Clusters of galaxies are formed in violent events, 
cluster mergers. These mergers can dramatically 
affect the observed properties of clusters, partic- 
ularly in the X-ray band (e.g., Markevitch et al. 
2002). Signatures of mergers include "cold fronts" 
(Markevitch et al. 2000; Vikhlinin, Markevitch, & 
Murray 2001b), merger shocks (e.g., Markevitch 



et al. 2002), and multiple cool cores (Markevitch, 
Sarazin, & Vikhlinin 1999). The application of 
simple hydrodynamical analysis to these merger 
signatures allows the merger Mach number and 
other aspects of the merger kinematics to be deter- 
mined (Vikhlinin et al. 2001b; Kempner, Sarazin, 
& flicker 2002; Sarazin 2002). In general, merger 
Mach numbers are modest, A4 ~ 2. 

Mergers also provide interesting tests of the role 
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of various physics processes in clusters. For exam- 
ple, very strong limits have been placed on the 
rate of thermal conduction and diffusion (Ettori 
& Fabian 2000; Vikhlinin, Markevitch, & Mur- 
ray 2001a; Markevitch et al. 2003b) and Kclvin- 
Hclmholtz instabilities and viscosity (Vikhlinin 
et al. 2001a). They also provide limits on the rel- 
ative roles of thermal versus nonthermal pressure 
sources in the intracluster medium (Markevitch 
et al. 1999). The role of merger shocks in accel- 
erating relativistic particles can also be studied 
(Markevitch & Vikhlinin 2001). Recently, the rel- 
ative distributions of gas, galaxies, and dark mat- 
ter in merging clusters have been used to limit 
the self-collisional cross-section of the dark matter 
(Markevitch et al. 2003a); collisional dark matter 
has been suggested to explain the mass profiles in 
dwarf galaxies and to resolve other problems (e.g., 
Spergel & Steinhardt 2000). 

On the other hand, merger effects can compli- 
cate the use of clusters of galaxies as cosmolog- 
ical probes. Clusters provide a number of im- 
portant cosmological tests, mainly because they 
arc the only objects which are large enough to 
represent a fair sample of materials in the Uni- 
verse, yet small enough to be relaxed (e.g., Bah- 
call 2000). Most analyses of clusters for cosmol- 
ogy assume that they are equilibrium systems; for 
example, the masses of clusters are often deter- 
mined from hydrostatic equilibrium, or more sim- 
ply by assuming an equilibrium relation between 
the cluster X-ray temperature or luminosity and 
the cluster mass. During major mergers, clusters 
are not in hydrostatic equilibrium. Also, compres- 
sion and heating during a merger can produce a 
large transient increase in the X-ray temperature 
and luminosity (flicker & Sarazin 2001). If equi- 
librium relations are applied to merging clusters, 
their masses can be greatly over-estimated (Ran- 
dall, Sarazin, & Ricker 2002). This can be par- 
ticularly important for high redshift clusters; in 
flux limited samples, merging clusters can be over 
represented due to their boosted X-ray luminosi- 
ties. It is important to understand these effects by 
studying relatively nearby clusters where more de- 
tailed dynamical information can be determined. 

fn this paper, we present new XMM-Newton 
observations of the Abell 1644 cluster. Our analy- 
sis of the Einstein observation suggested that this 
cluster was undergoing a major merger. Jones & 



Forman (1999) classified A1644 as a double clus- 
ter with roughly equal components based on the 
Einstein image. From a deprojection analysis of 
the Einstein data, White, Jones, & Forman (1997) 
found that A1644 harbors a small cooling flow of 
about 12 solar masses per year. Girardi et al. 
(1997) noted the possible presence of two com- 
ponents in the galaxy velocity distribution in the 
cluster, and listed two subclusters as part of this 
cluster. They also found that the cD galaxy had a 
large peculiar velocity, which suggests the possible 
presence of minor substructure. However, in the 
end they classified Abell 1644 as a regular clus- 
ter. No evidence for substructure has been found 
in recent optical and near-infrared observations by 
Tustin et al. (2001). A1644 is located close to a 
region of high galaxy cluster density. It is about 
15 degrees north (in Galactic coordinates) of the 
core of the Shapley supercluster (Fig. 1) at about 
the same redshift, z(A1644) = 0.0474 (Zabludoff 
et al. 1993). 




Fig. 1. — Location of A1644 relative to the core 
of the Shapley supercluster in Galactic coordi- 
nates. Shown are virial radii of a few clusters from 
Rciprich & Bohringer (2002) on top of a distribu- 
tion of 2MASS galaxies. 

2. Data Reduction and Analysis 

XMM-Newton, with its combination of high 
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throughput, broad spectral band, good spectral 
resolution, and fairly large field of view is ideal 
for studying this merging system. A 1644 was 
observed by XMM-Newton for a total duration 
of 23 ks during revolution 199 (observation ID: 
0010420201) on January 8, 2001. The EPIC-MOS 
cameras were operated in Full Frame (FF) mode 
and the EPIC-pn camera in Extended Full Frame 
(EFF) mode. The Thinl filter was used in all 
EPIC observations for protection against optical 
light contamination. Calibrated event files have 
been created with SAS V5.3. The calibration in- 
dex file has been created on May 7, 2002 based 
on the most up to date calibration files available 
at that time. Most further reduction and analy- 
sis steps have been performed using SAS V5.3 and 
the package HEAsoft V5.2. 

2.1. Flare Rejection 

XMM-Newton observations in general suffer 
from periods of very high background similar to 
Chandra observations (e.g., Lumb et al. 2002). 
There are several possible ways to 'clean' a data 
set of these flares. Our method is based on the 
fact that we will use independent blank-sky obser- 
vations (§ 2.3) to determine the main non-flare 
hard background component, the 'Particle In- 
duced Background' (PIB). (The PIB is distinct 
from the flares, which are believed to be pro- 
duced by 'soft' protons.) Thus, we must be able 
to clean the blank-sky background observations 
in the same way as the source observation. When 
determining the lightcurvc this requires ensuring a 
negligible 'contamination' by emission from astro- 
physical sources because otherwise the sensitivity 
to flare detection varies between source and back- 
ground observation. Here we choose for the MOS 
cameras the energy band 10-12 keV and for pn the 
band 12-14 keV similar to previous works (e.g., 
Majerowicz, Neumann, & Reiprich 2002), which 
are dominated by PIB (background induced di- 
rectly or through fluorescent instrumental lines by 
high energy particles) and soft protons (if present). 
Note that this approach has the drawback of low 
statistics and furthermore assumes that the flare 
characteristics in this energy band are representa- 
tive of those at lower energies as well. In § 2.3 we 
test the effect of increased statistics and the latter 
assumption. 

Times were binned in 100 s intervals, which ap- 



pears to be a reasonable compromise between good 
temporal resolution and low noise level. This time 
interval also is large compared to the read out cy- 
cles of the cameras. Further only events with pat- 
tern < 12 and flag = have been used to ensure a 
high rejection of potentially distorting events such 
as caused by hot/flickering pixels/columns or fea- 
tures close to CCD edges especially at the read 
out nodes. 

Inspection of the resulting lightcurves showed 
no obvious flares in the pn data (Fig. 2). MOS 
count rates were higher by up to a factor <~ 4 
in the beginning of the MOS observations com- 
pared to the remaining times (Fig. 3). MOS expo- 
sures started earlier than those for the pn camera; 
the lack of obvious flares in pn is consistent with 
the MOS lightcurves for the overlapping time in- 
terval. The obviously contaminated times in the 
MOS observations were removed by taking into ac- 
count only events with arrival times > 95360200 s 
(XMM-Newton time) for all further steps. 




9.5362e+07 9.5364e+07 9.5366e+07 9.5368e+07 9.537e+07 9.5372e+07 9.5374e+07 
TIME Is] 



Fig. 2. — Lightcurve for pn with 100 s binning. 

Before calculating mean count rates in order to 
test for background flares, conservative 'generic' 
cuts have been applied on the count rates, Cx- For 
pn (MOS) only times with 2 (1) < C x < 45 (25) 
counts per 100 s have been used. The lower thresh- 
old is especially important for pn data in order 
not to bias the calculation of the mean count rate 
by inclusion of counting mode times, times during 
which transmission of information for individual 
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Fig. 3. — Lightcurve for M0S1 with 100 s binning. 

events of one or more (most likely all four) quad- 
rants to the ground station does not occur. 

For the final cuts, upper and lower thresholds 
are required in order to treat statistical fluctua- 
tions in both directions symmetrically. Here we 
choose to define the good time intervals (GTIs) 
by including only times where the count rates lie 
within i 2d of the mean count rates. For the in- 
dividual detectors the count rate limits for GTIs 
are found to be 8-25 (pn), 3-15 (MOS1), and 2-15 
(MOS2) counts per 100 s. The effective exposure 
time (live time) after this procedure results in 11 ks 
for pn and 13 ks for MOS1 and MOS2. 

2.2. Vignetting Correction 

Depending on the nature of a detected event 
one should take into account vignetting caused 
by the baffles, mirrors, Reflection Grating Array 
(for the MOS cameras only), etc. in the analy- 
sis. For instance the X-ray photons from the tar- 
get source (A1644) and the Cosmic X-ray Back- 
ground (CXB) arc vignetted in the same way. The 
soft protons show a different vignetting character- 
istic (but hopefully we got rid of most of those in 
§ 2.1) and the PIB is not vignetted at all if the 
spatial variation of the fluorescent lines caused by 
the camera material itself is not interpreted as a 
vignetting characteristic. The noise in the CCDs, 
most easily visible as low energy events close to 
the CCD edges, is not vignetted, either. 



In general, however, it's not easily possible to 
determine the nature of an individual event. The 
combination of methods for flare rejection, vi- 
gnetting correction and background subtraction 
employed here aims at a proper treatment of all 
components for the spectral analysis. This strat- 
egy involves correcting all events for the vignetting 
expected for X-ray photons regardless of their na- 
ture at this point. 

The vignetting depends on detector position 
and energy. The vignetting correction has been 
accomplished with the SAS task evigweight which 
adds a weight column to the events file. The value 
of the weight is determined for each event by tak- 
ing into account the detector position and energy 
of the event. Weights are assigned relative to the 
on-axis effective area so the weights have values 
> 1. To accumulate images and spectra these 
weights are summed up, and errors are propagated 
including the vignetting correction. 

2.3. Background Subtraction 

In general, proper background subtraction is vi- 
tal for an accurate spectral analysis of extended 
sources. In our case, the cluster emission fills most 
of the fields-of-view of the detectors, so it is dif- 
ficult to determine the background from our ob- 
servations. Thus, we determine the background 
from so called blank-sky observations. The blank- 
sky observations are mostly high galactic lati- 
tude observations of regions with very few bright 
sources. Several such observations have been com- 
bined for the EPIC cameras mostly for the Thin 
filter in (mostly) FF mode (Lumb et al. 2002) 
to achieve good statistics and reduce the effects 
of cosmic/sample variance. We have made sure 
that for the energy range and signal-to-noise ra- 
tio in spectral fits here, the differences in Charge 
Transfer Inefficiency (CTI) corrections between 
the pn FF and EFF mode are not significant. The 
blank-sky observations are available at the XMM- 
Newton calibration web page 1 . 

The background observations have been recast 
on the sky in order to match the source obser- 
vations using the SAS task attcalc 2 . This allows 
using sky coordinates for regional selections and at 

1 http: / / xmm.vilspa.esa.es / external/xmm_sw_cal/calib/epic_filcs.shtml 

2 See Andy Read's scripts at 
http://www.sr.bham.ac.uk/xmm3/ . 
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the same time estimating the background from the 
same detector regions as the source, which in turn 
ensures that the significant variation of the fluores- 
cent line emission with detector position is prop- 
erly taken into account. By treating the blank-sky 
backgrounds identically to the data, the effect of 
incorrectly applying a vignetting correction to the 
PIB (in the source and background observations) 
is removed. 

The same flare rejection method (§ 2.1) has 
been applied to the background observations. The 
count rate limits for GTIs here are found to be 8- 
31 (pn) and 2-16 (MOS1 and MOS2) counts per 
100s. 

The PIB shows variation of about 8% (la) as 
well as a long term decrease of about 20% between 
the beginning of XMM-Newton observations and 
January 2001 (Katayama et al. 2002). In order 
to account for these variations, total field-of-view 
count rates of source and background observations 
have been calculated using the same energy, pat- 
tern, and flag selection criteria as for the GTI de- 
termination (§ 2.1). The source to background 
count rate ratios 0.99 ± 0.03 (MOS1), 0.92 ± 0.03 
(MOS2), and 0.88 ± 0.02 (pn) have been found 
(1(7 errors), consistent with the expected varia- 
tion. For the calculation of these PIB scaling 
factors the events have not been corrected for vi- 
gnetting. 

As mentioned in § 2.1, cleaning of flares in the 
high energy bands generally used does not neces- 
sarily imply that all flares have been removed. Ac- 
tually, we did find some weak remnant flares in the 
blank-sky observations when plotting lightcurves 
(GTI screened as described in § 2.1) in broad en- 
ergy bands (see also Lumb et al. 2002). We there- 
fore tested the variation in (10-12 kcV) source 
to background count rate ratios when calculating 
GTIs in the energy bands 7-12 and 0.3-12 keV for 
MOS1 and MOS2. We found that the variation is 
less than the (1 a) statistical uncertainty and con- 
clude that the influence of remnant flares is negli- 
gible for calculation of background scaling factors 
here. 

To check whether the limitation to high en- 
ergy bands for the calculation of background scal- 
ing factors affects the results significantly we have 
calculated them also using the energy ranges 0.3- 
14kcV (pn) and 0.3-12 keV (MOS). In this case, 
of course, only the CCD parts shielded from X-ray 



photons (and soft protons) are usable. We found 
results consistent with the total field-of-view re- 
sults for pn and MO SI and a slightly larger value 
for MOS2 (1.00 ± 0.02). In order to exclude the 
possibility that an imperfect correlation between 
the normalization of the continuum and the flu- 
orescence lines (dependent on detector position) 
for the PIB caused this slightly larger than ex- 
pected variation, we also used the energy range 
2-7 keV to calculate the background scaling fac- 
tors. The results were consistent with the broader 
energy range used previously. 

In summary, following these tests the back- 
ground scaling factors 0.96 (MOS) and 0.88 (pn) 
have been assigned to the observations. Most of 
the tests above yield results consistent with the 
expected statistical uncertainties. To account for 
the difference between MOS and pn scaling factors 
a conservative 8% systematic error has been added 
to the full energy range of all background spectra. 
Note also that it is implicitly assumed here (and in 
most other works that follow similar procedures) 
that the spectrum of the PIB does not vary when 
its normalization changes. Currently, there seems 
to be no strong evidence for significant spectral 
variability (e.g., Katayama et al. 2002) apart from 
the fluorescence lines. 

The main drawback of this method is that in 
general one gets the background contribution from 
the CXB wrong. The PIB dominates the hard part 
of the total background but the CXB is impor- 
tant for the lower energies. Ideally one would esti- 
mate the difference in CXB contribution to source 
and background observations from a region with 
negligible source emission in the source observa- 
tion or another observation very close by (e.g., 
Markevitch & Vikhlinin 2001; Pratt, Arnaud, & 
Aghanim 2001; Majerowicz et al. 2002; Arnaud 
et al. 2002). This is not feasible in the case of 
A 1644 since cluster emission fills the entire field 
of view and the only other XMM-Newton observa- 
tion within 5 degrees to date is likewise filled en- 
tirely with NGC 5044's emission. Here, we model 
a possible CXB difference in the spectral fitting 
procedure. 

In order to estimate the spectral shape of a pos- 
sible CXB difference between the source and back- 
ground spectra we have extracted spectra from 
an irregular region with decreased cluster emis- 
sion in the north-western part of the observation 
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and fitted a two temperature model. A high tem- 
perature component which is supposed to account 
for the cluster emission and a low temperature, 
zero redshift, solar abundance component which 
should represent the CXB excess relative to the 
background observations. The model fits indicate 
that in this region the cluster emission still dom- 
inates. The ratio of the number of cluster pho- 
tons to CXB photons equals about 4 in the en- 
ergy range 0.3-10 keV and about 2 in the energy 
range 0.3-0.5 keV. A best fit value - 0.2 keV is 
found for the low temperature component. This 
estimate is consistent with values from studies of 
soft excess emission in other clusters (e.g., Kaas- 
tra et al. 2003). The origin of this component is 
uncertain; Kaastra et al. (2003) argue that it is 
due to intcrgalactic gas located around clusters. 
Our aim here is not to try to discuss the source of 
this soft CXB excess, but rather to determine its 
contribution to the cluster emission and account 
for it in the model fits (§ 2.4). 

2.4. Spectral Fitting 

Spectra from MOS1, MOS2, pn have been 
fit simultaneously. The same pattern and flag 
selections have been applied as for the genera- 
tion of GTIs, except that only single and double 
events have been used for the pn camera. The 
response matrices were mll_r7_im_alL2000-H- 
09.rmf, m21_r7_im_all_2000-ll-09.rmf, and epn_cf2 
0_sdY9.rmf for MOS1, MOS2, and pn respectively 
in addition to on-axis ancillary response files cre- 
ated with the SAS task arfgen. 

A comparison of individual pn single and dou- 
ble event spectral fits gives consistent results. 
The response of the pn camera depends on de- 
tector position. However, it has been verified 
that for the accuracy needed here, use of the 
epn_ef20_sdY0.rmf response matrix does not result 
in significant changes of best fit parameter values. 
Individual fit results from the three cameras are 
consistent with each other. Point sources and hot 
pixels have been excluded. 

XMM-Newton's point spread function (PSF) is 
much improved compared to the PSFs of previ- 
ous satellites capable of spatially resolved spec- 
troscopy, e.g., ASCA. Nevertheless, for objects 
with strong temperature and surface brightness 
gradients its effect can become important. How- 
ever, the main emphasis in this work is on the 



comparison of different segments of spectra at the 
same radii from cluster centers. In this case radial 
PSF effects should not cause any artificial differ- 
ences. 

Spectral fits have been performed within XSPEC 
using a model (wabs*mekal)+mckal. Parameters 
for the second thermal component have been fixed 
at a redshift of zero, a metal abundance of A = 1 
solar (Anders & Grevesse 1989), and tempera- 
ture of kTx = 0.2 keV, as discussed in § 2.3. The 
normalization is determined from the fits. 

Source spectra have been regrouped to have 
at least 50 counts per bin. Spectral fits have 
been performed nominally in the energy range 0.3- 
10.0 keV for MOS and 0.3-12.0 keV for pn. All er- 
rors are quoted at the 90% confidence level for one 
interesting parameter unless noted otherwise. 

3. Results 
3.1. Imaging 

Figure 4 shows the combined MOSl-MOS2-pn 
image of nearly the full XMM-Newton field-of-view 
of A1644. The cluster has a very complicated sur- 
face brightness distribution on all scales. A main 
cluster to the southwest and a smaller sub cluster 
to the northeast are easily identified. The emis- 
sion surrounding both core regions is highly non- 
spherical. There is stronger emission to the south 
of the sub clump than to the north. The core re- 
gion of the main clump itself contains a displaced 
core-within-a-core. 

The Digitized Sky Survey (DSS) image in Fig. 5 
shows that the X-ray peaks of main cluster and sub 
cluster coincide with two Bright Cluster Galax- 
ies (BCGs). The relative line-of-sight velocity 
component of the two central galaxies is about 
600 km s -1 . The X-ray contours to the north of 
the sub clump appear to be slightly compressed. 
Radio observations of the two BCGs show that 
they are both radio sources (e.g., Owen & Ledlow 
1997). 

The regions around the two BCGs are shown 
in more detail in Figs. 6 and 7. Included are X- 
ray contours of a 19 ks Chandra observation taken 
on the same day as the XMM-Newton observa- 
tion and available from the public archive. These 
contours illustrate the complexity of the emission 
in the very cluster cores. Note that also in the 
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Fig. 4. — Background, exposure, and vignetting 
corrected, adaptively smoothed, combined MOS1- 
MOS2-pn count rate image of A1644 in the energy 
range 0.3-5.0 keV (R.A. and Dec. in equatorial co- 
ordinates, epoch J2000.0). The linear artifacts are 
due to an imperfect exposure correction near the 
chip boundaries. 



Fig. 6. — The DSS image of the inner region of 
the main clump, overlaid with the X-ray contours 
from Chandra (black) and, in the very center, the 
1.4 GHz radio data (red; from M. Ledlow). 




J 94.900 



Fig. 5. — X-ray contours from the image in Fig. 4 
overlaid onto an optical DSS image. 



Fig. 7. — The DSS image of the inner region of 
the northeast sub clump overlaid with the X-ray 
contours from Chandra (black) and, in the very 
center, the 1.4 GHz radio data (red; from M. Led- 
low). 



very center of the sub clump the X-ray contours 
to the north appear slightly compressed (Fig. 7). 
This apparent compression is on a much smaller 
scale than the compression already noted in Fig. 5. 
Optical, radio, and X-ray centers coincide almost 
exactly for both sub clusters. 

3.2. Spectroscopy 

The intraclustcr gas temperature is a funda- 
mental observable for a gravitational cluster mass 
determination based on the hydrostatic assump- 
tion. The overall best fit temperature based 
on a large elliptical region centered at R.A. = 
194.356 deg and Dec. = — 17.363deg encompass- 
ing both sub clusters (semi major axis = 733", 
semi minor axis = 528", position angle = 45 deg) 
is kT x = 3.83 ± 0.06 kcV. The best fit metal abun- 
dance is 0.32±g;g| solar for the same region. 

Note that the column density of neutral Galac- 
tic hydrogen, iVjj, has been fixed at the value 
inferred from 21 cm radio measurements for our 
Galaxy (5.33 x 10 20 cm~ 2 ; Dickey & Lockman 
1990) for all results. In general this value is 
marginally larger but consistent with best fit val- 
ues from direct spectral fits. Best fit values for 
temperatures and abundances change only well 
within their uncertainties if the column density is 
treated as a free parameter. 3 

Figures 8 and 9 show the radial temperature 
profiles for the main and sub clump, respectively. 
Note that the center for the innermost bin of the 
main clump is offset from the center of the other 
bins in order to account for the core-within-a-core 
structure. The reduced \ 2 values (and the best fit 
iVn values if left as a free parameter) for the inner- 
most bin of both sub clusters are too large (Tab. 1 
and 2). This is likely caused by higher tempera- 
ture gas in the line of sight. For these two regions 
we therefore included a third thermal component, 
i.e., the model (wabs*(mekal+mckal))+mekal was 
fitted. The limited statistics require freezing the 
temperatures and abundances of the third compo- 
nents to the values found further out. We therefore 
fixed the abundances at 0.5 solar and the tempera- 



3 Note that neglecting the CXB difference between source 
and background observations, i.e., fitting a model 
wabs*mekal to the spectra here, results in values sig- 
nificantly lower than the Galactic value, or, if Nn is fixed 
to the Galactic value, to artificially decreased temperatures 
and increased \ 2 - 



tures at 3.5 and 3 keV for the main and sub clump, 
respectively. As a result the best fit central tem- 
peratures decrease further for both main and sub 
clump (open diamonds). 
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Fig. 8. — Radial temperature profile for the main 
clump. The open diamond shows the best fit tem- 
perature if a third thermal component is included 
(see text). 



4 r 



1 — o - 
^ 2: 



100 200 300 

Radius [arcsec] 

Fig. 9. — Radial temperature profile for the sub 
clump. The open diamond shows the best fit tem- 
perature if a third thermal component is included. 

Each of the two temperature profiles appears 
surprisingly similar to temperature profiles of re- 
laxed, apparently undisturbed clusters: a drop in 
the center to about 1/3 of the ambient gas temper- 
ature, an isothermal structure in the outer parts, 
and weak indications for a slight temperature drop 
in the very outermost regions accessible. 

Note that the addition of a third component is 
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only a crude approximation to account for effects 
caused by projection, PSF (see, e.g., Markevitch 
2002; Majerowicz et al. 2002), and a possible real 
multiphase temperature structure. The reduced 
X 2 values of the central bins become significantly 
smaller but are still large even after addition of 
the third component. Also, these effects need not 
be limited to the central bins alone. A1644 is an 
extremely irregular object (see Figs. 4 and 12). 
This fact makes A1644 so interesting but at the 
same time complicates deprojection attempts, for 
which spherical symmetry usually is assumed. In 
this work, we therefore focus on projected tem- 
perature estimates and qualitative comparisons to 
temperature maps of simulated clusters. 
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Fig. 10. — Metal abundance profile for the main 
clump. The open diamond shows the best fit abun- 
dance if a third thermal component is included. 

The corresponding metal abundance profiles 
are shown in Figs. 10 and 11. In general, the 
abundances increase towards the centers of the 
two clumps, although there seems to be a drop 
in abundance in the innermost region of the sub 
clump. However, this abundance drop vanishes 
when performing the fit with three thermal com- 
ponents (open diamond). This drop can therefore 
be regarded as artificial, even though it should be 
noted that the precise values depend on the tem- 
perature and abundance values adopted for the 
third component. 

It is worth noting that weak indications for 
higher normalizations of the CXB excess emission 
towards the center of the main clump are found 
which may support the findings of Kaastra et al. 
(2003), although a discussion of the nature of the 



1.0 : 
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1 0.6^ 

< 0.4 ; 
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Radius [arcsec] 

Fig. 11. — Metal abundance profiles for the sub 
clump. The open diamond shows the best fit abun- 
dance if a third thermal component is included. 



CXB excess is beyond the scope of this paper (see 
§ 2.3). 

Figure 12 shows a hardness ratio (HR) map. 
To create this map, first, source photon images 
from MOS1, MOS2, and pn (using only single 
events) in the energy band 2.0-10.0 keV have been 
created and combined. Then, the correspond- 
ing exposure maps have been created and com- 
bined, which also account for vignetting. Using 
the combined image and exposure map, an adap- 
tively smoothed and exposure and vignetting cor- 
rected source image has been created with the 
SAS task asmooth (signal-to-noise ratio = 15). 
Subsequently, a similarly created (though simply 
smoothed with a = 7.5), exposure and vignetting 
corrected, combined background image has been 
subtracted from the source image. The same pro- 
cedure has been followed to create an image in the 
0.3-2.0 keV energy band, using the same smooth- 
ing templates as created for the 2.0-10.0 keV im- 
age. The ratio of the latter two images is shown in 
Fig. 12. Soft emission is shown in dark and hard 
emission in bright. We have verified that using 
the source smoothing template also for smoothing 
the background image does not result in significant 
changes in the relevant region. Note that the CXB 
excess emission has not been subtracted in the HR 
map. This may bias the HR map in the outermost 
regions of the observation. However, we will find 
below that spectral fits that take the CXB excess 
into account correspond well to the features seen 
in the HR map. Therefore, it is not expected that 
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Table 1 

Fit results for the main cluster profile 



Radial range ["] 


fcTx [kcV] 


A 


red. x 2 


D.O.F. 


0-30 
0-30 a 
30-77 
77-201 
201-325 
325-449 


2 44+ - 16 

i qc+0.15 
1 " 3o -0.13 
o 91 +0.12 

4 2S+ 013 

4 39+ 15 
^• oa -0.15 

4 15+0.19 


0.85+° ;?° 
94+ 1 ' 11 

-0.36 
n 7 o+0.09 

n on+0.05 
u -°'-'-0.05 

24+ - 05 
u - z4 -0.05 

21+ 06 
u - zi -0.06 


1.80 
1.46 
1.19 
1.11 
1.09 
0.97 


103 
103 
312 
769 
808 
872 



a Third thermal component applied (see text). 



Table 2 

Fit results for the sub cluster profile 



Radial range ["] 


fcT x [kcV] 


A 


red. x 2 


D.O.F. 


0-30 


, 9 o+0.06 
± - zo -0.06 


n 9 o+0.07 
U -^°-0.05 


2.65 


54 


0-30 a 




55+ - 50 
u - oo -0 .29 


1.75 


54 


30-77 


r, oc+0.21 

z,oo -0.19 


44+ ' 14 
u -^-0.12 


1.27 


108 


77-201 


o 7n +0.19 
a - ' U -0.17 


49+ - 1 " 
u - y,a -0.09 


0.96 


404 


201-325 


3.98t°;i 


0.29±0;°| 


0.89 


570 



a Third thermal component applied (see text). 
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the region shown is biased significantly. 




Fig. 12. — Hardness ratio map of A1644, based on 
the ratio of the combined MOSl-MOS2-pn count 
rates in the energy bands 0.3-2.0 keV and 2.0- 
10.0 keV. The map has been corrected for back- 
ground, exposure, and vignetting, and is adap- 
tively smoothed. Soft emission appears dark and 
hard emission bright. Also shown are surface 
brightness contours from Fig. 5 and regions se- 
lected for spectral analysis. 




OTn ian £□ 1a. □ 1s:5rna.a 

Fig. 13. — Temperature map of A1644. Also 
shown are surface brightness contours from Fig. 5. 

In Fig. 13 an adaptively binned temperature 
map is shown. The map has been created using 
all three detectors and making use of the Interac- 
tive Spectral Interpretation System (ISIS; Houck 
& Denicola 2000). The pixel size is 14.4". Pixels 
are nearly independent at the center of the main 
clump but are correlated in other regions. The 
actual spectral extraction region used for the fits 
depends on the surface brightness; the largest ex- 
traction region is l'on a side. Most pixels shown 
black had fewer than 800 source counts using the 
maximum smoothing scale. A few more pixels 
with negative error > 30% were masked out. The 
abundance of the main mekal component has been 
fixed here at the average value 0.32 solar (§ 3.2). 
The features seen in the HR map (to be discussed 
in detail below) seem to correspond to features in 
the temperature map. This indicates the consis- 
tency of both approaches and we'll use the HR 
map, which is more closely related to the data, as 
guidance below. 

The hardness ratio map and the X-ray im- 
age suggest a number of interesting regions in 
the cluster which may have differing temperatures 
and/or abundances. In general, clusters as ir- 
regular as A 1644 often show more complex tem- 
perature structure. Therefore, we also extracted 
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spectra in segments of annuli for the main clump. 
The radial boundaries for the nine segments of the 
main cluster correspond to the ranges for the three 
outermost annuli given in Tab. 1. Figure 12 shows 
these and other selected regions as well as surface 
brightness contours overlaid onto the HR map. 
(Note that the cluster emission extends beyond 
the outermost shown surface brightness contour.) 

Best fit temperature and metal abundance val- 
ues for the segments are shown in Figs. 14 and 
15 (and Tab. 3), respectively. It is apparent that 
there is significant nonradial temperature struc- 
ture, and that the relatively regular radial tem- 
perature profile in Figure 8 is the result of averag- 
ing over annuli with large azimuthal temperature 
variations. In particular, the temperatures of the 
segments to the west of the main clump (R3, R6, 
and R9) are all significantly lower than all regions 
to the east (Rl, R2, R4, R5, and R7) except for 
R8. These findings correspond well to those in the 
HR map (Fig. 12). 



5.5 




3.0 L 1 
2 4 6 8 10 

Region 



0.6 
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2 4 6 8 10 
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Fig. 15. — Metal abundance determinations of 
lected regions of the main clump. 
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Fig. 16. — Gas temperature determinations of 
lected regions of the sub clump. 



Fig. 14. — Gas temperature determinations of se- 
lected regions of the main clump (see Fig. 12). 

We also determined temperatures and abun- 
dances for more complicated regions selected as 
being of interest based on the XMM-Newton im- 
age and the HR map. This allows us to search for 
cool/metal-rich trails (dark in the HR map) or hot 
spots (bright in the HR map) . Whether brightness 
fluctuations in the HR map correspond to signifi- 
cant temperature/abundance changes is evaluated 
by direct spectral fits of selected regions. Note 
that not all artifacts, e.g., inexact exposure cor- 
rection close to CCD chip boundaries, have been 
removed in the HR map in Fig. 12. The spectral 



0.8 - 
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Fig. 17. — Metal abundance determinations of 
lected regions of the sub clump. 
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analysis, however, is not affected by this. 

The region to the south of the sub clump (R3) 
appears fairly soft, whereas a small region to the 
southeast (R4) appears hard. The spectral fit re- 
sults (Fig. 16 and Tab. 4) reveal that indeed the 
temperature determination for R3 is significantly 
lower than the value for the rest of this annulus 
(R2, which excludes Rl) and especially than that 
for R4. Furthermore, the metal abundance of R3 
appears enhanced (Fig. 17). The best fit tem- 
perature for R5 (the small dark tail to the south 
of R3 in Fig. 12) is similar to that for R3. On 
the other hand, Rl, the region north-northeast of 
the sub clump, has a high temperature similar to 
R4. 

4. Discussion 

We first consider some implications of these ob- 
servations for the use of clusters of galaxies as cos- 
mological probes. How much do the X-ray flux es- 
timates based on resolved observations differ from 
simple estimates based on global measurements, 
such as those for clusters only observed with the 
ROSAT All-Sky Survey (RASS)? How much do 
mass estimates differ when only a broad beam 
overall gas temperature estimate is available? 

The overall flux from A1644, which extends 
beyond XMM-Newton's field-of-view, is /x(0.1 — 
2.4 kcV) = 4.03 x KT 11 ergs s" 1 cm~ 2 (Reiprich & 
Bohringer 2002) based on RASS data. (A1644 was 
not observed in pointed mode with the PSPC.) It 
is worth noting that A1644 has been treated as 
a single cluster in the HIFLUGCS (Reiprich & 
Bohringer 2002) and REFLEX (Bohringer ct al. 
2001) catalogs, primarily due to the limited ex- 
posure time (< 300 s) of the RASS at A1644's 
position and the relatively small projected sepa- 
ration of the two cluster components. The XMM- 
Newton image indicates that the X-ray flux (fx) 
ratio between the main and sub clump is about 
3:1. This means instead of one cluster with fx ~ 
4 x 10~ n ergss -1 cm~ 2 one actually has two clus- 
ters with about fx ~ 3 and 1 x 10~ n ergs s _1 cm~ 2 
each. This could affect luminosity and mass func- 
tions significantly if such major mergers were com- 
mon. The effect of an incorrect treatment would 
be to artificially increase the number of high mass 
clusters and decrease the number of low mass clus- 
ters. Using the shape of the mass function to con- 



strain the cosmic; mass density. f2 m , and the am- 
plitude of mass fluctuations, expressed as erg, indi- 
vidually would then result in artificially low values 
for f2 m and high values for a%. With the ongoing 
analysis of the complete sample of HIFLUGCS 
clusters with Chandra and XMM-Newton, we will 
be able to quantify the fraction of such clusters 
and their influence on derived cosmological pa- 
rameters, as well as to study the difficulty of es- 
timating survey selection functions in the vicinity 
of bright clusters (not only due to double clusters 
but also line of sight projections) in general. 

A simple broad beam temperature estimate 
may be biased low compared to the ambient gas 
temperature due to cool emission in the dense 
(high emissivity) cores. For instance, the temper- 
ature in a large annulus around the main clump 
in A1644 (201-325") is 4.40 keV. If we take this 
as the ambient temperature which is characteris- 
tic of the cluster gravitational potential, it is a 
factor of 1.15 higher than a broad beam temper- 
ature estimate including both clumps (3.83 keV). 
Since M cx T 1 - 5-2 - (e.g., Finogucnov, Reiprich, & 
Bohringer 2001), this factor translates into an un- 
derestimate of a factor of 1.23-1.32 in the cluster 
mass. 

The surface brightness structure of both sub 
clumps is obviously messy. However, the radial 
temperature profiles (Figs. 8 and 9) are smooth 
and appear rather similar to profiles of relaxed 
clusters (e.g., compare to Fig. 1 in Allen, Schmidt, 
& Fabian 2001). This might suggest that A1644 
is moderately relaxed, and that the temperature 
structure may not be strongly affected by the in- 
teraction of the two sub clumps. However, a 
more detailed examination of the thermal struc- 
ture shows that the cluster is not at all relaxed. As 
the HR map (Fig. 12) indicates and Figs. 14 and 
16 confirm, the temperature structure is also quite 
complex. The region between the two clumps is 
significantly hotter than the region to the west of 
the main clump (the region which appears least 
disturbed in X-ray surface brightness). This may 
indicate that the gas in the region between the 
clumps has been heated up by adiabatic compres- 
sion or shocks. Also, the core-within-a-core struc- 
ture of the main clump (Fig. 4) may be caused by 
core oscillations induced by the passage of the sub 
clump (e.g., Tittley & Henriksen 2003) or possibly 
by a remnant of a previous merger. 
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Table 3 

Fit results for the main cluster regions 



Region feT x [keV] A red. x 2 D.O.F. 



1 


4 71+0.28 
4 - ' ^o^ 


32+° 10 
U - JZ -0.10 


1.06 


259 


2 


4 54+ ' 27 

^■°^-0.26 


Q no + 0.10 


1.08 


267 


3 


q q 7 +0.19 


41 +°- 08 


0.96 


439 


4 


4 70+ ' 27 
4 -' u -0.24 


n 99 +0.09 
u - zz _0.08 


1.07 


289 


5 


4 89+ - 33 


17 +0.09 


0.94 


316 


6 


o O4+0-27 


27+ - 10 

u - z( -0.09 


1.11 


410 


7 


4 69+ - 25 


n 99 +0.08 
u - z -0.08 


0.83 


318 


8 


4 11+0.27 
4 - 11 -0.27 


0.24+^° 


1.03 


340 


9 


3 64+ 33 
C> - D -0.31 


o.is±g:ii 


0.93 


413 



Table 4 

Fit results for the sub cluster regions 



Region feT x [kcV] A red. x 2 D.O.F. 



1 


r 24+1-79 


7n +1 ' 34 

u - ' u -0.70 


0.89 


30 


2 


n 7 o + 0.31 
d -'°-0.29 


36+ - 13 


0.92 


217 


3 


3 18+ - 24 
■J- 1 °_0.22 


"59+ - 14 
u - oa -0.13 


1.27 


165 


4 


5 49+ 1 ' 45 

J-i»_0.99 


42+ 69 
u -* z -0.42 


1.46 


31 


5 


3.09l°;« 


31 +°- 24 
U - O1 -0.22 


1.28 


45 
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Furthermore, both the surface brightness and 
temperature structure around the sub clump are 
peculiar. Based on the X-ray image, the gas ap- 
pears to be slightly compressed to the north and 
clearly elongated to the south and west. This sug- 
gests that this sub clump is moving through the in- 
tracluster medium (ICM) of the main clump, and 
losing some fraction of its gas in a trail behind 
the core of the sub clump. Given this dynamical 
picture what temperature structure might be ex- 
pected? One may expect heated gas in front of the 
sub clump (projected roughly to the north) due to 
adiabatic compression or shocks. One might also 
expect the trail to the south to consist of some- 
what cooler, more metal-rich gas stripped from 
the sub clump. This appears to be consistent both 
with the HR map and with the extracted spectra 
in the sub clump regions Rl, R3, and R5. 

Figs. 5 and 7 show that the X-ray peak of the 
sub clump and the optical center of the bright- 
est galaxy coincide. This suggests that the sub 
clump has not passed through the core of the main 
clump; if it had, the very high ram pressure would 
have stripped the sub clump core, which would lag 
behind the collisionless dark matter and galaxies 
(e.g., Markevitch et al. 2002). The angle of the 
cool, metal-rich gas trail also suggests an off-center 
collision. One concern with this picture is that the 
dense gas in the core of the sub clump would then 
be interacting with lower density gas in the outer 
parts of the main clump. Would one expect the 
dense sub cluster to lose an observable amount of 
gas during its travel through the ICM? 

We have compared A1644 to a hydrodynamical 
simulation of a small cluster undergoing an off- 
axis merger with a larger cluster. This simulation 
was done completely independently; no effort was 
made to adjust the initial conditions to match the 
observation of A1644. The temperature map of 
this simulation at a time when the smaller clus- 
ter has passed by the core of the main cluster for 
the first time is shown in Fig. 18. Comparison to 
the HR map of A1644 (Fig. 12) shows an interest- 
ing similarity. There is hotter gas in front of the 
smaller cluster and cooler, ram-pressure stripped 
gas behind. This trail of stripped gas cools adia- 
batically as it expands. There is a hotter region to 
the southwest, which is similar to region R4 seen 
in A 1644. The larger scale Tx distribution around 
the main cluster in the simulation is also similar 



to that observed in A1644. There appears to be 
slightly cooler gas to the southeast in the simula- 
tion due to stripped gas from the smaller cluster, 
which agrees with the R8 region in the HR map 
and spectra. 4 




Fig. 18. — Mass weighted temperature map from 
a hydrodynamical simulation. Temperature in- 
creases from black to red to yellow. Also shown 
are X-ray surface brightness contours. 

The indication for an increased metal abun- 
dance in the cool trail of A1644's sub clump makes 
other scenarios like cooling of intracluster gas onto 
the moving sub clump (as seen in A1795 on much 
smaller scales, Fabian et al. 2001) more difficult to 
reconcile with the data. Note that the radio source 
in the center of the sub clump (Fig. 7) may pro- 
vide additional energy for gas removal from central 
denser regions. 

5. Summary 

XMM-Newton imaging spectroscopy of the 

4 An animation of this simulation is available at 
http://A1644.dark-energy.net. Note that cool 

trails from sub clusters are commonly seen in 
simulations. Compare, e.g., to the simulations 

by Motl et al. (2003, see also the movie at 
http: / /www. astro. Virginia. edu/coolflow/abs.php?regID=184). 
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galaxy cluster A1644 has shown two sub clus- 
ters undergoing an off-axis merger with remark- 
able detail. The X-ray analysis is also supported 
by optical data on the cluster, and the resulting 
dynamical picture agrees very well with an in- 
dependent hydrodynamical simulation of such an 
off-axis merger. The findings imply that we see 
gas that has been removed from the smaller sub 
cluster possibly by ram pressure from the motion 
through intracluster gas of the main clump, pos- 
sibly augmented by the effect of the central radio 
source of this sub clump. 

The code to create the X-ray temperature map 
made use of the Interactive Spectral Interpretation 
System (ISIS; Houck & Dcnicola 2000). M. Led- 
low made available the radio images. This work 
was supported by NASA XMM-Newton Grants 
NAG5-10075, NAG5-13088, and NAG5-13737. T. 
H. R. acknowledges support by the Celerity Foun- 
dation through a Post-Doctoral Fellowship. The 
XMM-Newton project is an ESA Science Mis- 
sion with instruments and contributions directly 
funded by ESA Member States and the USA 
(NASA). The Digitized Sky Surveys were pro- 
duced at the Space Telescope Science Institute un- 
der U.S. Government grant NAG W-2166. The 
images of these surveys are based on photographic 
data obtained using the Oschin Schmidt Telescope 
on Palomar Mountain and the UK Schmidt Tele- 
scope. This publication makes use of data prod- 
ucts from 2MASS, which is a joint project of the 
University of Massachusetts and IPAC/Caltech, 
funded by NASA and NSF. 
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